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Abstract 

We find the finite volume QCD partition function for different quark masses. This 
is a generalization of a result obtained by Leutwyler and Smilga for equal quark masses. 
Our result is derived in the sector of zero topological charge using a generalization of the 
Itzykson-Zuber integral appropriate for arbitrary complex matrices. We present a conjec- 
ture regarding the result for arbitrary topological charge which reproduces the Leutwyler- 
Smilga result in the limit of equal quark masses. We derive a formula of the Itzykson-Zuber 
type for arbitrary rectangular complex matrices, extending the result of Guhr and Wettig 
obtained for square matrices. 
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1 Introduction 



In general, numerical simulations offer our only access to the QCD partition function. For 
sufficiently small volumes, however, the partition function is dominated by the constant 
modes which makes analytic treatment possible. (See, e.g., ref. [Jj].) In a similar spirit, 
Leutwyler and Smilga identified a parameter range within which the mass dependence 
of the partition function is completely determined by the underlying structure of broken 
chiral symmetry. This range is given as 

t « L « — , (1) 
A m % 

where L is the linear size of the 4-dimensional Euclidean box, A is a typical hadronic 



mass scale and m n is the mass of the Goldstone modes (i.e., ~ ymK for quark mass 
m). The lower limit of this range ensures that the partition function is dominated by the 
Goldstone modes. The upper limit ensures that these modes can be treated as constant 
modes. If chiral symmetry is broken according to SU(Nf) x SU(Nf) — > SU(Nf), the 
QCD partition function in the range ([]]) and for vacuum angle 9 is given by [[J 

Z(M,6)= ( dUexp(vZRe(e i9/N nrUM)) . (2) 

JU£SU{N f ) V ' 

Here, S is the chiral condensate, and M is the mass matrix which can be taken as diagonal 
without loss of generality. In the sector with topological charge v, the partition function 
is given as 

Z„{M) = I** dde~ ive Z(Q) . (3) 

Z7T JO 

For equal quark masses, Leutwyler and Smilga obtained an exact analytic expression 
for this partition function 

Z v {m) = det I v+j ^{mV E) , (4) 

where % and j run from 1, • • • , Nf. The importance of this partition function lies in the 
fact that it enables us to determine the volume dependence of the chiral condensate. This 
is of interest in lattice QCD simulations, where the volume is necessarily finite. 

Recently, lattice QCD calculations have been performed to determine the connected 
and disconnected contributions to the scalar susceptibility |§. The determination of 



these susceptibilities requires differentiation of the partition function with respect to two 
different masses. In order to determine volume dependence of the scalar susceptibility in 
the range (|l|), we require a generalization of (§) to different quark masses. This is the 
primary objective of the present paper. 

In section 2, we analyze the two flavor case and generalize the result to an arbitrary 
number of flavors. In section 3, we present a generalization of the Itzykson-Zuber formula 
valid for arbitrary rectangular complex matrices. For square complex matrices, our result 
reduces to that obtained by Wettig and Guhr [|]. In section 4, we derive the finite volume 
partition function for different masses in the sector of zero topological charge. We also 
make a conjecture of the result for arbitrary topological charge and show that it reduces 
to ((H) for the special case of equal quark masses. 



2 The finite volume partition function for Nf = 2 and 
its generalization to arbitrary Nf. 

For two flavors and in the range (|I|), the QCD-partition function is known for different 
quark masses. For vacuum angle 9, it is given as 

2(fl) = vrA(vs/i) , (5) 

where I\ is a modified Bessel function and where 

/i 2 = m\ + m\ + 2mim 2 cos(6>) . (6) 

The partition function in the sector with topological charge v can be obtained by inte- 
grating over 9 according to (H). Remarkably, this integral can be expressed analytically. 
After some manipulations the result can be written as 

W = 2) = - J ^- J det 

where 



I v {xi) 



X2K 



(7) 



x k = m k VJ:. 
3 



(8) 



This suggests the generalization to three flavors 

I v { Xl ) x 1 I' v {x 1 ) x\I'l{x x ) 

* (9) 



16 

Zy{N f = 3) = 2 ,w 2 2W T 2 _ ~2\ ^ et 

l x 2 x lA x 2 x 3Jl x 3 x lJ 



4(a; 2 ) x 2 I' v {x 2 ) xll'l(x 2 ) 
I v (x s ) x 3 I' v (x 3 ) xll>l{xz) 



where the numerical prefactor is chosen such that for equal masses the result of Leutwyler 
is reproduced. The generalization to an arbitrary number of flavors is now obvious. We 
define a Vandermonde determinant as 

A(x 2 ) = U(xl-xf), (10) 

k<l 

and an Nf x Nf matrix as 

A M = x^I^ixk) k,l = l,---,N f . (11) 
The partition function is then given as 



7 < \ n detA 

1 1 A(x z ) 



(12) 



where the normalization constant 

C Nf =2 N ^ N f~ i y 2 Y[(k-l)\ (13) 



Nr 



k=l 

is determined by the limit of equal quark masses. 

In section 4, we shall prove this formula for the special case v = 0. It will be shown 
that this formula reduces to the Leutwyler-Smilga finite volume partition function for 
arbitrary v. 

3 The Itzykson-Zuber integral for complex rectangu- 
lar matrices 

In this section we offer a derivation of the extension of the Itzykson-Zuber integral to the 
case of arbitrary complex rectangular matrices using the diffusion equation method. The 
result for square matrices has also been obtained in 0. Our derivation is patterned on 
the argument for Hermitean matrices, which has been discussed widely in the literature. 
(See, e.g., f| 0, f|.) The expression which will be required in section 4 for the calculation 
of the finite volume partition function for different masses is given in (|33|). 



Let a and p be arbitrary complex (i.e., non-Hermitian) Ni x N2 matrices. Without 
loss of generality, we assume that u = Ni — N 2 > 0. Any such matrix can be diagonalized 
in the form |§ 

a = U^SV (14) 

where 

«-(o)- < 15 > 

and where the square diagonal matrix 5 = diag(si, • • • , sjva) nas nonnegative real entries. 
The matrices (U,V) parameterize the coset space U(Ni) x U(N2)/[U(1)] N2 . Similarly, p 
can be written as p = U'^RV. 

In the following, we evaluate the integral 

(7r J Nl N 2) J MU, V) exp (-±tr[(<7 - p)\a - p)}) (16) 

using the diffusion equation method and exploiting the invariance of the measure, dp(U, V), 
which is taken to be the Haar measure of U(Ni) x U(N2)/[U(1)] N2 . The function 

Hp, t) = ( nt y NlNa) J M {-fA^ - p) V - p)]) <p(°-) (17) 

with integration measure d[o\ = llm=i Hn=i dRea mn dIma mn satisfies the diffusion equa- 
tion 

Ni N 2 q2 a 

EEsir'te') = s «. us) 

m =l n =l u h'mn^Hmn ul 

As initial condition, we choose an invariant function, i.e., f{<j) is a function of the eigen- 
values of cr only. Then, using the invariance of the measure, we find immediately that 
F(p, t) is a symmetric function of the eigenvalues of p. 

In order to proceed, we express both the Laplacian and the integration measure in the 
'polar coordinates' introduced by the diagonalization (|i~4|), 



d[a] = Qd[S] j 2 {S) dp(U, V) (19) 

where d[S] = lln=i ds n . Here, the constant f2, which depends on the convention adopted 
for the measure of the group, will be fixed later. The Jacobian is given by J(S) = j 2 (S) 



with 

N 2 

m = n 8m- n wa{&) , (20) 

n=l 

and A(S 2 ) is the Vandermonde determinant defined in (|T0"|). Because F is an invari- 
ant function, it satisfies a diffusion equation which involves only the radial part of the 
Laplacian: 

S^^V^ar • < 21 > 

This equation can be simplified materially with the introduction of a reduced wave func- 
tion, 

f(R,t) = j(R)F(R,t) . (22) 



Now, (|2lT) reduces to 



N 2 / Q2 



<9 2 

dr 2 " 7 



f{R,t) = 4?l . (23) 



Remarkably, this differential equation is separable.^] Performing the differentiations of 
j{R), it can be rewritten as 

Ni ( d 2 4z/ 2 - 1 1 \ df 

Because of the presence of the factor j(R) in (|22|), f(R,t) is an antisymmetric function 
of the eigenvalues. Therefore, the solution of (gj) is given by an integral over a Slater 
determinant 

f(R,t) = Jd[S\^d^\g(r k ,ar,t)\ j(S)<p(S) , (25) 

where g(r, s; £) is the kernel of 

d 2 4z/ 2 - 1 dg 

d^ 9 ~^^ 9 = ^ • (26) 
The kernel can be expressed in terms of the regular eigenf unctions of the Bessel equation 
of order v: 



u 



" + [k 2 - (4z/ 2 - l)/4r 2 ]u = . (27) 



1 In (|], a separable equation was obtained for iVi = N2 by the substitution of A(R 2 )F(R 7 t) instead 
oij(R)F(R,t). 



This is obtained following a separation of variables which leads to a time dependence of 
the form exp(— k 2 t/4). The regular eigenfunctions are given by 

u fc (r) ~ Vkr J v {kr) , (28) 

where J v is a Bessel function. Recalling the orthogonality relation for Bessel functions, 
the kernel of (p6) can be written as 

g(r,s;t) = 6(t) dkk^Tse- kt/A JJkr)JJks) . (29) 
Jo 

This can be evaluated to give 

^(r.ajt) = -v^exp f j — \l u (—J (30) 

where is a modified modified Bessel function. 

Finally, we equate the definition of F(p, t) in ( |TTD with its expression in terms of the 
kernel (|30|) as given by ( p2|) and (^5|). Since this equality is valid for an arbitrary choice 
of the initial condition <p(a), the integrands of d[S] must be the same. 

Hence, we arrive at 



J dfi(U, V) exp (-\tr[(a - p)\a - p)]) 



t N 1 N 2 ~N 2 2 N2 7T NlN2 1 

H e t 



exp — — n-H 



/A(^)A(i? 2 ) . (31) 



« = T^lv^T- 7T7 - (32) 



Here, the value of f2 follows from the normalization integral calculated in pO 

n- v iii!0' + ^- 1)' 

where we have used the convention that J dp(U, V) = 1. 

This result enables us to calculate the Itzykson-Zuber integral for arbitrary complex 
matrices 

/ d/*(IW ^p (RetrCrtSra) = C Nl C Na , 2 ,^ det | 1^)1 /A(S 2 )A(R 2 ). 

(33) 

When z/ = 0, the product of factorials in this expression is understood to be 1. The 
constants Cn 1 and Cn 2 can be evaluated to be 

n 

C n = 2™(™- 1 )/ 2 J[(k- !) ! • (34) 
fe=i 

In the special case iVi = N2 for which i/ = 0, our expression reduces to the result of Guhr 
and Wettig || apart from a normalization constant. 



4 The finite volume partition function for different 

masses 



The finite volume effective partition function of QCD in the sector with topological charge 
v, defined in and fl2|), can be written as an integral over U(Nf) instead of SU(Nf) 

Z u = ( dn(U)(detUy exp(vi: Re tr MU^) . (35) 

JU(Nf) V ' 

For v = 0, this result can be obtained from (|33| ) by taking R = VY>M and S = 1^. 
Because of the singularity as S — > 1, the final result requires a careful analysis of this 
limit. 

In (|33"D, we take S = 1 + 5S and expand the modified Bessel functions to order (5S) Nf . 
The result is 

N f j'-i 

WW) = £ Jjh^ I^Xn) Ssl 1 + 0(5S N f) , (36) 
which can be written as the product of the matrix A defined in (|TTD and the matrix 

The determinant of B can be written as 

det(fl) = -A(5S) = 77— A ((1 + 5S) 2 ) , (38) 



with the normalization constant defined in (0). Hence 

det A 

which is simply the result conjectured in (O) above. This result is now proved for v = 0. 

We have not proved the result ([12]) for arbitrary v. However, we offer one non-trivial 
check of this conjecture by demonstrating that (|IIJ) reduces to the finite volume partition 
function of Leutwyler and Smilga in the limit of equal masses. We start with the expression 

A k3 = C l ^rl l I^\r k ). (40) 

Using the recursion relation 

^ = W + ~h (41) 
dr r 



and adding a suitable multiple of the column to the left of the column in question, we 
arrive at the matrix 



A, - ^ r l,Nl r j ~ l T 



N f 'k 



Jk) 



(42) 



This is evidently correct in going from the first to the second column and, hence, true in 
general. The fact that the coefficient of I u+ \ in (f41~|) is 1 guarantees that the determinant 
will not be affected by this rearrangement. 

To realize the limit r = mVT,, we write = r + 5rk and expand each element 

in A in a Taylor series through order 5r^ f 1 . The result is that 



Akj — C Nf 1 ^2 



r j 1 h +j -i{r)\ 5r p k 1 



(p - 1)! drP- 1 

This can be written as the product of the matrix with elements 

dP- 1 



(43) 



M. 



jp 



C 



A' 



(44) 



and the matrix Bj P = (Sr^Y 1 /{p — !)!• As in (^3), we have 



N, 



det(B) = J] 



-A(8R) 



i n\N f (N f -i)/2 



A((R + 5R)' 



The determinant of the matrix M can be simplified by using the recursion relation 

dl v _ v 

—r~ — lv-i J-u 

dr r 



(45) 



(46) 



and adding a suitable multiple of the row immediately above the row in question. As 
before, this rearrangement does not alter the determinant. This leaves us with 



M, • -> r 1/Nf r i- l T ■ , 

lvl kj ' ^Nt ' 1 v+]-k ■ 



(47) 



Now, the factor r 7 ' can be extracted from each column and used to eliminate the r depen- 
dence in the prefactor. Thus, we arrive at the final result 



det(A) = det(M) det(S) = det !„+,•_*. 

k,j 



(48 



which is precisely the result in || as given in 



5 Conclusions 



We have obtained the finite volume QCD partition function for different quark masses in 
the range 1/A C L C l/m n . This result generalizes the finite volume partition function 
obtained previously by Leutwyler and Smilga for the case of equal quark masses. In order 
to derive this result, we were led to generalize the Itzykson-Zuber integral to arbitrary 
rectangular complex matrices. The integral for square matrices, first obtained by Guhr 
and Wettig, leads immediately to the proof in the sector of zero topological charge. Based 
on the result for two flavors and the general result for v = 0, we have conjectured the 
result for arbitrary v and Nf. As a decidedly nontrivial check of this conjecture, we have 
shown that the result of Leutwyler and Smilga emerges in the limit of equal quark masses. 

We wish to note a remarkable coincidence. Consider the Itzykson-Zuber formula for 
complex rectangular matrices (with v equal to the difference between the number of rows 
and columns) in the same limit considered for v = 0. Up to a factor, this leads to the 
finite-volume partition function in the sector of topological charge v. While we can offer 
no explanation of this coincidence, it may be useful to note that, in the chiral limit, the 
joint eigenvalue density of the random matrix model associated with the finite volume 
partition function depends only on the combination v + Nf [ llj • Clearly, more work lies 
ahead in this area. 

Acknowledgements 

The reported work was partially supported by the US DOE grant DE-FG-88ER40388. 



References 

[1] P. van Baal, Nucl. Phys. Proc. Suppl. 20 (1991) 3. 
[2] H. Leutwyler and A. Smilga, Phys. Rev. D46 (1992) 5607. 
[3] F. Karsch and E. Laermann, Phys. Rev. D50 (1994) 6954. 
[4] C. Itzykson and J.B. Zuber, J. Math. Phys. 21 (1980) 411. 

10 



[5] T. Guhr and T. Wettig, preprint |hep-th/9605110 . 
[6] M.L. Mehta, Comm. Math. Phys. 79 (1981) 327. 
[7] T. Guhr, J. Math. Phys. 32 (1991) 336. 

[8] P. Di Francesco, P. Ginsparg, and J. Zinn- Justin, Phys. Rep. 254 (1995) 1. 

[9] L. Hua, Harmonic Analysis of Functions of Several Complex Variables in the Classical 
Domain, trans. L. Ebner and A. Korani, American Mathematical Society, Providence, 
RI, 1963. 

[10] T. Nagao and M. Wadati, J. Phys. Soc. Japan 60 (1991) 3298. 

[11] J.J.M. Verbaarschot, Phys. Rev. Lett. 72 (1994) 2531; Phys. Lett. B329 (1994) 351; 
Nucl.Phys. B427 (1994) 434. 



11 



